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I. INTRODUCTION 

Even though quantum chromodynamics (QCD) over 
time has become the generally accepted theory of the 
strong interaction, it has stiU proven difficuh to use this 
theory to derive low-energy hadronic observables such as 
masses and decay constants. Lattice QCD simulations 
are an alternative approach where the functional inte- 
grals are evaluated on a discretized spacetime through 
numerical Monte Carlo simulations. At present however, 
computational limitations has hindered such simulations 
for light particles, due to the fact that they can prop- 
agate over large distances, requiring very large lattice 
sizes. Because of this, most simulations have been per- 
formed with heavier quark masses than those of the phys- 
ical world. Typically, the quark masses used in present 
lattice simulations fulfill mu,d > rris/S. In order to get 
physically relevant predictions, these results then have to 
be extrapolated down to the physical masses of ~ ms/25. 

Chiral perturbation theory (xPT) Q, |M |1| is an effec- 
tive field theory which provides a theoretically correct 
description of the low-energy properties of QCD. If the 
quark masses of the lattice simulations are low enough 
so that the corresponding xPT calculations can be con- 
sidered accurate enough, this allows for a determination 
of the low-energy constants (LECs) of xPT by means of 
a fit to the lattice simulations. This in turn makes es- 
timates of hadronic low-energy observables possible. In 
particular, masses and decay constants of the physical 
pseudoscalar mesons can be determined in this way. 

It has however proven difficult to reach the chiral 
regime, and therefore many lattice simulations have been 
performed with so-called quenched QCD. In this the- 
ory, the effects of closed (sea) quark loops have been 
neglected, since such loop effects require repeated evalu- 
ation of fermion determinants, which is computationally 
extremely expensive. Another option is to not neglect 
these loop effects altogether, but instead to introduce a 
separate quark mass mgoa for the sea-quarks loops. Such 
procedures are referred to as partial quenching (PQ), and 
lead into a space of unphysical theories. This has the 
advantage over full QCD calculations that results with 



more values of the valence quark masses can be obtained 
with a smaller number of values of sea quark masses since 
varying the latter is computationally more expensive. 

Since unquenched QCD may be recovered from par- 
tially quenched QCD (PQQCD) in the limit of equal 
sea and valence quark masses, it follows that QCD and 
PQQCD are continuously connected by the variation of 
sea-quark masses. In contrast, this is not true for fully 
quenched theories. 

It has been shown that xPT can be extended to include 
both quenching and partial quenching 0, 0, 0, Q . For 
partial quenching this is particularly interesting, because 
it allows for determination of the physically relevant 
LECs of xPT by fits of partially quenched xPT (PQxPT) 
to partially quenched lattice simulations (PQQCD), see 
e.g. the discussion in Q. The reason for this is that the 
formulation of PQxPT is such that the dependence on 
the quark masses is explicit, and thus the limit of equal 
sea and valence quark masses can also be considered for 
PQxPT. It follows that xPT is recovered as a continuous 
limit of PQxPT, just as QCD is from PQQCD. In partic- 
ular, the LECs of xPT, which are of physical significance, 
can be obtained directly from those of PQxPT. Clear dis- 
cussions of this point as well as calculations at next-to- 
leading order (NLO) and references to earlier work can be 
found in the papers of Sharpe and Shoresh |^ . There 
are two variants of PQxPT possible, with and without 
the supersinglet $0 degree of freedom. This degree of 
freedom is the partially quenched analog of the singlet rj' 
and is expected to be heavy compared to the other pseu- 
doscalars for not too large quark masses. We work in this 
paper with the variant without the $0 ■&] and with three 
quark flavors. 

In previous papers we calculated the masses and de- 
cay constants of the charged, or off-diagonal in fla- 
vor, pseudoscalar mesons in PQxPT to next-to-next- 
to- leading order(NNLO) or 0{p^) in the chiral count- 
ing 0, ^> ^1 ■ For these mesons and their masses 
the needed constants of PQxPT to NNLO could be de- 
termined by varying quark masses in the valence sector. 
This way the structure of the extrapolation to the phys- 
ical limit of equal valence and sea quark masses and low 



2 



values for the light quark masses can be studied in detail. 

Here we now turn to the neutral meson sector. The 
structure of the neutral or flavor-diagonal propagator be- 
yond lowest order (LO) in PQ^PT is discussed in [^l^. 
It is shown there that the fully resummed case has double 
poles at the masses related to the valence sector as well 
as single poles at masses related to the valence sector and 
to the neutral meson masses in the sea quark sector. The 
neutral meson masses in the valence sector are the same 
as the equal quark mass limit of the charged valence sec- 
tor meson masses, i.e. the mass in the valence qq channel 
is the same as in the valence qq' channel with niq = ruqi . 
This is certainly not true in full QCD. A consequence is 
that the neutral meson masses for full QCD are not a 
simple limit of the neutral meson masses in the valence 
sector of PQQCD but one has to study directly the me- 
son masses in the sea-quark sector. At NLO this can be 
seen by the fact the LEG Lj cannot be determined from 
the masses in the valence sector but it is known to play a 
role in the eta mass Q. Even in the sea quark sector, 
can only be determined from the neutral meson masses 
if at least two different sea quark masses are used, i.e. 
m ^ rris, since the contribution of Ly to the eta mass is 
proportional to (m — mg)^. This means that determining 
the chiral behavior and checking general convergence for 
neutral meson masses is quite difficult in lattice calcula- 
tions. 

Sharpe and Shoresh showed in Sect. VI in 8] that at 
NLO in PQxPT a measurement of the ratio of correlators 

r, ^ {^r^^it,p^ 0)T^jj{x ^ 0)) , . 

°^ (7r,,(t,p = 0)7r,,(x = 0)) 

allows to determine L7. Tr^j is the source for a pseudo- 
scalar meson with quark flavors qiq^ and we take the va- 
lence quark masses rrii = rrij but i and j different flavors. 
For large times t this correlator behaves as 

Vt 

i?o(t-^oo) = ^, (2) 

with Mij the mass of the charged meson with quantum 
numbers qiqj. The quantity V is thus a physical quan- 
tity measurable in PQQCD on the lattice. The large t 
behavior in Q follows from the double pole structure 
of the full neutral propagator . In |^ |3| 2? was also 
studied to NLO in PQxPT. There they showed that 
can be determined from P, even in the case with all sea 
quark masses equal. This way the LECs relevant for 
the eta mass can be determined from PQ^PT, including 
the more detailed studies allowed by varying the valence 
quark mass. 

V vanishes when the valence and the sea quark masses 
are equal. The main part of this paper is devoted to 
studying T), or more generally, the full neutral propaga- 
tor, to NNLO in PQxPT. We therefore first show how the 
double pole structure from the full neutral meson propa- 
gator follows from an all order resummation of diagrams 



rather than the arguments used in [3, |9j . We have per- 
formed this resummation for the general mass case. In 
the main text we present the result only for the simplest 
case. The most general case can be found in App. IbI 
Afterwards we calculate all relevant parts to NNLO and, 
in particular, we show the results for V. We also discuss 
how V allows to determine all needed LECs to NNLO for 
the eta mass. 

This paper is organized in the following manner: Sec- 
tionlTTlintroduces the PQxPT formalism and an overview 
of the notation used for loop integrals and quark masses. 
Section IIIII contains a derivation of the resummed neu- 
tral propagators and the double-pole coefficient V. Sec- 
tion ^3 describes how to calculate V to NNLO and also 
contains the analytical expressions. Section presents 
a numerical analysis of the results. An extension of the 
fitting strategies for extracting the various LECs from 
lattice QCD calculations is then found in section IVII 
followed by a brief discussion of our conclusions in Sec- 
tion |VlIl 

II. xPT AND PQxPT 

In this section we give a very short overview of some 
aspects of xPT and PQxPT u p to NNLO. Lectures on 
standard xPT can be found in [ijj NNLO results there 
are reviewed in 0. A good discussion of PQxPT can 
be found in |^ and the NNLO aspects are discussed in 
our earlier papers, but most extensively in We refer 
to those references for more details. 

In PQxPT, a mechanism which gives different masses 
to sea quarks and valence quarks is introduced in a sys- 
tematic way by adding to xPT explicit sea quarks, as 
well as unphysical bosonic ghost quarks. The latter can- 
cel exactly all effects of closed valence quarks due to their 
different statistics, provided that their masses are iden- 
tical to those of the valence quarks. This results in a 
modification of the chiral symmetry group of xPT. In 
PQxPT, the symmetry group is given by 

G = SU (n^al + nsce,\n^al)L X SU{n^i^l + risoa|"-val)_R, • (3) 

The precise structure of G is somewhat different as dis- 
cussed in 1^ . This theory contains rival valence and 
rival ghost quarks, as well as rigca flavors of sea quarks. 
As the PQ theories include bosonic ghost quarks, they 
are not normal relativistic quantum field theories since 
they violate the spin-statistics theorem. However, under 
the assumption that the low-energy structure of such a 
theory can be described similarly to the case of normal 
QCD, one arrives at an effective low-energy theory in 
terms of a matrix U, according to 

U = exp(iV2^/F^ . (4) 

The matrix $ contains the Goldstone boson fields gen- 
erated by the spontaneous breakdown of the chiral sym- 
metry group Q to its diagonal subgroup. $ has a more 
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complicated flavor structure than in ordinary xPT be- 
cause of the different types of quarks present. In terms 
of a sub-matrix notation for the flavor structure 



UaUb Uadp UaSb 
QaQb = I daUb dadp daSb 
SaUb Sadb SaSb 



(5) 



where we have used three quark flavors u, d and s in each 
sector, the matrix $ becomes 



qvqv 



qsqv 



\ qsqv 



qvqs 



qsqs 



qsqs 



qvqs 



qsqs 



qsqB 



\ 



(6) 



where the labels V, S and B stand for valence, sea and 
bosonic quarks, respectively. The size of each sub-matrix 
depends on the exact number of quark flavors used. 

The quarks qv, qs and their respective antiquarks are 
fermions, while the quarks qs and their antiquarks are 
bosons. Thus a combination of a ferniionic quark and a 
bosonic quark yield a fermionic (anticommuting) field, 
while a combination of two fermionic or two bosonic 
quarks result in a bosonic field. Each sub-matrix in 
Eq. © therefore consists of either fermionic or bosonic 
fields only. The bosonic quarks are given the same masses 
as the corresponding valence quarks in order to cancel 
the contributions from closed valence quark loops. The 
above formalism is often referred to as supersymmetric 
PQxPT. 

The Lagrangian structure of PQxPT is the same as for 
n-flavor xPT, provided that the traces of matrix products 
in those Lagrangians are replaced by supertraces. The 
supertraces are defined in terms of ordinary traces by 



Str 



A B 
C D 



TrA-TrD, 



(7) 



where A, B,C and D denote block matrices. For exam- 
ple, the block B corresponds to the [qvqB] and [qsqs] 
sectors of the field matrix in Eq. © and contains anti- 
commuting fields, while the block D represents the [qBqs] 
sector of Eq. © . A more detailed discussion about the 
Lagrangians and LECs can be found i n 1131 . The ex- 
pressions for the Lagrangians are in 0, Il7|. The LECs 
at NLO are labeled L[ and those at NNLO . The 
subtraction scale dependence has been suppressed in this 
notation. 

The version of PQxPT considered in this paper has 
three flavors of valence quarks, three flavors of sea quarks, 
and three flavors of bosonic 'ghost' quarks. The different 
quark masses are identified in the following calculations 
by the flavor indices i = 1,...,9, rather than by the 
indices u,d,s and V,S,B of Eqs. ((SJl and ©. The re- 
sults are expressed in terms of the quark masses via 
the quantities Xi = "^Boiriqi such that Xi;X2,X3; belong 



to the valence sector, X4iX5iX6 to the sea sector, and 
XTjXSiXg to the ghost sector. The latter ones do not 
appear in the results since the ghost quark masses are 
always set equal to the masses of the corresponding va- 
lence quarks, such that X7 = Xi,X8 = Xi and X9 = Xz- 
Xi is the lowest order meson mass squared for a charged 
pseudoscalar meson with two quarks with masses rrii. 

For the explicit NNLO calculations in this paper we 
only consider the case with all valence quark masses 
equal, xi = X2 = X3j and all sea quark masses equal, 
X4 = X5 = X6- Thus only the two quark-mass param- 
eters xi and X4 appear in the final results, but we also 
introduce the combinations 



Xl4 = (xi +X4)/2, ^u = Xi-X4. 



(8) 



Xi4 is the lowest order meson mass squared for a pseu- 
doscalar meson with two quarks with masses mi and 771,4. 
Rf^ is the only one of the many combinations of quark 
masses needed in the more general case Jl^, J,3] . 

The analytical expressions depend on several one- and 
two-loop integrals. After carrying out the regularization 
and renormalization, finite contributions remain, which 
are defined as 

^(x) = -71-16 xlog(x/M^), 

Xi log(Xi/M^) - X] log(Xj7M^) 



(5(X,X,X;0) = -7ri6/(2x). 



X* - X3 



(9) 



where /i denotes the renormalization scale and Trig = 
l/(167r^). For the case of equal quark masses the ex- 
pression for B reduces to 

B(X,X;0) = -Trie (1 + log(x/M')) . (10) 

Some combinations of these functions arc naturally gen- 
erated by the dimensional regularization procedure. We 
therefore define 



A{x\e) 
B{x,x;0,e) 



A{x)V{2nie x) + ^16 X (^Vl2 + 1/2), 
A{X)B{X, X; 0)/(^i6 x)-^(x)V(27ri6 x') 
+Tri6 (^Vl2 + 1/2). (11) 



Evaluation of the two-loop sunset diagram introduces an- 
other class of integrals, whose finite contributions have 
been denoted by iJ^, ijf , H^^ and ijf , H^^ . The 
latter three functions are derived from the former three 
by differentiation with respect to before going on-shell. 
These functions all contain an integer argument which is 
used to distinguish between the 8 possible configurations 
of double and single poles in the integrands. A detailed 
discussion can be found in for PQxPT, and in p^ . 
where the single-pole version of these integrals first ap- 
pear in the context of ordinary xPT. 
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III. THE RESUMMATION 
A. Introduction 

The resummation procedure is analogous to the re- 
summation of the charged mesons [lO, 11, 12, 13J, except 
that the neutral sector has mixing. One solution to this 
problem is to first diagonalize the lowest order and then 
do the resummation that way. That was done for the 
isospin breaking in the 7r° and rj masses in [T9l |. Doing 
the resummation that way requires that the lowest or- 
der propagator matrix in the neutral sector is invertible. 
In usual xPT this is done by working in the lowest or- 
der eigenstate basis, i.e. one has a two by two matrix 
structure in the neutral sector for three quark flavors. 

In the PQ theory, the same mechanism could be used, 
but this is complicated by the fact that one has many 
more states and the double poles in the lowest order prop- 
agator as well. It is therefore easier to work in the flavor 
basis rather than the lowest order eigenstate basis. Af- 
ter the removal of the $0 degree of freedom this means 
that the lowest order propagator matrix is not invertible. 
We thus have to perform the resummation in a different 
way. Sharpe and Shoresh ^ did this by keeping the 
$0 degree of freedom and only removing it after the re- 
summation. Here we remove it from the start and show 
that we reach the same conclusions. 

We start by writing the general propagator as 

G,,ki{x) = {n\T {^,jix)^kim \n) , = qrQj, (12) 

where Q denotes the vacuum of the interacting the- 
ory. Similarly, we denote the sum of all one-particle- 
irreducible (IPI) diagrams with incoming quantum num- 
bers as $y and outgoing as by —iT^iju- Note that 
this notation differs slightly from that of the previous ar- 
ticles 0,0,0,0]. For the neutral sector we define the 
9x9 matrices in the {VV , SS, BB) basis as the subsets 
of Gijki and Y^i^a with 



Gujj{x)^{n\T{^.,,{x)^,,m\n). 



(13) 



Except for Add. IXI we will work with the Fourier trans- 
form of these quantities. 

We now look at all contributions to the full propaga- 
tors. First we recapitulate the charged propagator, i ^ j, 
case. There we have that 



G 



ijkl — Gijji 



(14) 



because of the conservation of flavor. The lowest order 



propagator G^, 



ijkl 



and the self-energy E. 



ijkl 



have 



the same flavor structure as Eq. (|14|) . The sum over 
all diagrams for the charged propagators as depicted in 
Fig. ^becomes a diagonal sum with no internal indices to 
be summed over. This resummation of the diagrams in 
terms of lowest order propagators and self-energy con- 
tributions becomes a simple geometric series. When 



summed, it gives a shift to the position of the pole of 
the lowest order propagator, see e.g. 0|. The position 
of this resummed pole defines the physical masses. 

For the neutral sector the procedure is similar, but 
less straightforward due to the mixing. Neither the self- 
energy nor the lowest order propagator is diagonal, we 
thus need to keep the nine by nine structure as defined in 
Eq. (|13|l . Consequently, the resummation has the struc- 
ture 



ij 



Gij + G%{—i)T,kiG'}j 
+Gik{-i)^kiG°„^{-i)T,jnnG°j 



, (15) 



as depicted in Fig. ^ If — lEG*' had been invertible it 
would have been a simple geometric series which could be 
resummed via matrix inversion as done e.g. in [T9l |. We 
perform below this resummation explicitly taking into 
account the exact structure of the matrices G° and — lE. 
The simplest case with all valence quark masses equal 
and all sea quark masses equal is treated here, and is re- 
ferred to as the (H-l) case following the notation used in 
our earlier work. The generalization to the more general 
cases can be found in App. 



B. Matrix Structure of G° and -iE. 

The lowest order neutral propagators have the general 
form nil 111 



G"(fc) = G'iJk)S,,-GUk)/r 



(16) 



The full structure can be found in Add. iBland {l3i |. 

For the (1-f 1) mass-case we are interested in, the ex- 
pression simplifies considerably. By taking the appropri- 
ate mass limits it becomes 

G° = il® H-iK(g)D, (17) 

where I is the 3x3 unit matrix, K is given by 

1 1 1 



K = i I 1 1 1 
3 1 1 1 1 



(18) 



and the two matrices H, D are given by 



H 



/ a 











/3 






u 





-a / 





a -f 7 a a + ■f 
a f3 a \ . (19) 

a + J a a + 7 



The symbols a, f3, 7 denote the different types of propa- 
gator terms which remain after the mass limit has been 
taken, and are explicitly 



1 



-Xi 



1 



Xl - Xi 



I and K give the structure inside the various sectors 
as the matrix shown in Eq. jSJ . H and D give the sector 



k ^^^V 1 j i k ^^^V 1 ni^^^V n j 

+ I I f I + 



m ^ ^ n ^\ P J 




+ • • • • 



FIG. 1: The diagrams contributing to the two-point function GJ^ . The Unes are lowest order propagators. The filled circles 
correspond to the sum of all one-particle-irreducible diagrams, the self-energy. The labels near the self-energy show the flavor 
index to be summed over. 



structure as in Eq. 0. The contributions with I corre- 
spond to the connected diagrams and those with K to the 
disconnected diagrams when thinking in terms of quark 
lines. Similar remarks apply to all the other matrices 
discussed below. 

Thus, one has a rather clean matrix structure for G^, 
involving just the three parameters a, /3,7. Furthermore, 
since K is a projection operator, = K, one avoids 
arbitrarily high powers of K. 

Next, we investigate the matrix structure of the self- 
energy terms — zE. In order to do this, we again turn 
the attention to the propagators Gijki- For the case 
of one flavor of valence quarks and one flavor of sea 
quarks, Sharpe and Shoresh showed that several 
constraints on the structure of Gijki follow from the (su- 
per)symmetries inherent in the PQ theory, and that the 
argument can be generalized to the more general cases. 
In App. we show this extension expicitly. From this 
derivation follows that the block structure of the neutral 
sector, before the $o degree of freedom has been inte- 
grated out from the theory, can be written as 



G = 




(21) 



with 



s = si , 
t = tK. 



r = rK 
u = s'l 



r'K. 



(22) 



The symbols s, r, t, s' , r' denote unknown functions, 
whose form we are not interested in for this discus- 
sion. We only need the structure of their arrangement 
in Eq. (|?T|l . Defining, 



W 
S 



T ■ 
R 



(23) 

s"i(ts"Hu"i -rs^i), (24) 



one then finds that the matrix 

S 



R 



G- 




(25) 



is an inverse to G. When checking this explicitly, it is 



useful to note that us tu = ts 



This matrix is a 



straightforward generalization of the inverse given for the 
1+1 flavor case in l^. 

If we now decompose G~^ as 



(26) 



where (G")^^ is the inverse of the lowest order propa- 
gator matrix, then — iS corresponds to the sum of the 
one-particle irreducible self-energy diagrams. Since G*^ 
has the same block structure as G, the inverse {G^)~^ 
will have the same block structure as G"^ as well. Thus, 
subtracting (G°)~^ from G~^, one finds that — iE has the 
same block structure as well, namely 



E 




(27) 



where R,S,T,W are different matrices from those in 
G~\ but with the same matrix structure. The same 
structure should prevail after removing the <i>o degree 
of freedom |^ . One can write E in a convenient di- 
rect product-notation as well. From the definitions of 
i?, 5", T, W one can rewrite them as 

R ^ pK, S^al, (28) 
T = tK, W ^wil + W2K. (29) 

where p,a,T,wi,W2 are as of yet unknown functions. 
With this notation, E can finally be written as 



1(g) A + K(g)B, 



(30) 
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where 



















Wl 


• 






^0 











(31) 



Similar arguments show that a is also the self-energy 
needed for the off-diagonal or charged meson with both 
quarks having the same mass. 



The resummation 



(2') 

Similarly, the r,\ ' combine to 



(3) 

and, since 7 has a double pole, the r„ sum up to 



.(3) 



-iiXi - Xi) 
In all, we thus have the result 



(39) 



(40) 



We now perform the resummation of Eq. (|15(l . Since K 
is a projection operator, a typical term in the sum looks 
like 

G"(-iEG'°)" = iG°(l(g)(AH)" + K(K){all other terms}) . 

(32) 

Both A and H are diagonal, the first part thus rcsums 
component by component as a geometric series, giving 



/ 



G 



AH 











V 











(33) 



The "K-part" seems less trivial. Writing out a few of 
the terms in a symbolic manipulations program such 
as Maple, one quickly sees a recurring structure which 
suggest a simple formula for the general term. Having 
this formula, a proof that it holds for all orders then 
follows by induction. Denoting each 3x3 block of G 
by = 1,2,3 we now list the resulting expres- 

sions for each sector and resum them. First, the sectors 
(1, 1), (1, 3), (3, 1) and (3, 3), i.e. r, all result in the same 
expression for a typical term giving 



na{p + Wl +W2 — 2r) — {n + l)a{a + 7) 



a a , 



where n = 0, 1, 2, 
into three parts 



For convenience, now divide r„ 



' n ' n ' n i 



(35) 



with 



.(1) ^ 

n 

.(2) _ 



.(3) _ 



-i7(n 4- l)(a-ck)" . 



For rn\ note that 



1 



ri=0 



(n-f-l)a;". 



(36) 



(37) 



Using this and the definition of a from Eq. H2U|) , it follows 



that the Vn'^ resum into 

r(i) = (p + wi+W2-2r) 



(P^ - Xi - cry 



(38) 



J5 - Xl - 

i{p + Wl + W2 - 2t - a - (xi - X4)) 

(p^ - xi - 



(41) 



The other sectors are not needed for this paper, since they 
do not contain any double pole contributions. However, 
we still list the full result here for completeness. Pro- 
ceeding in the same way, we get that the neutral sector 
resummed propagator becomes 

G = iI(g>lis-iK(g>'Ds, (42) 

where the two matrices Hs, Dg are given by 












• 













-a / 





a + J a a -f 7 

a P a 
a -f 7 a a + J 



(43) 

The symbols a, /?, 7 denote the various resummed prop- 
agator terms, given by 



7 



p2 - Xl - cr ' ^ - X4 - Wl 

Xi - X4 + o- + 2r - p - Wl - W2 
(p2 - xi - cr)^ 



(44) 



We see that the full propagator has exactly the same 
structure as the lowest order propagator in this simplest 
mass case. The resummed propagator matrix involves 
the three parameters a,/3,7, which are the direct coun- 
terparts to the lowest order parameters a,/3, 7 of G^. 



D. V from G 

We now expand the double pole around the position of 
the pole. The position is the same as the charged mass 
squared since a is also the self-energy for that case. The 
position of the zero of — xi — cr we thus refer to as 
M^i^. It is the physical mass of the charged mesons (in 
the degenerate mass-case, there is only one such mass). 
We expand around the double pole term and get 



-iZV 



(45) 
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The renormalization factor 2, is defined as 



= 1 



(46) 



Here and below we have indicated the dependence on the 
momentum and the Xi of the various quantities. They 
also depend on the LECs of PQxPT. The definition of V 
in 1)45(1 coincides with the one given in Eq. (0). 
Eqs. H42|) and (|45|l then imply that T) is given by 



V = 



Xi - X4 + o- + 2r - p - wi - W2 



(47) 



Note that the factor 1/3 from K has been included here 
as well. 



where a' contains derivatives of the self-energy dia- 
grams of O(p^), and cr'^^^ derivatives of those of 0{p^). 
Note that differentiation lowers the order of the diagram 
contributions by two, and therefore these terms are of 
O(p^) and 0{p^) respectively. It is sufficient to use the 
lowest order mass instead of M^^ in a'^^^ since the dif- 
ferentiated diagrams in that term are already of 0{p'^). 
However, in the case of ct'^'*-' the argument M^^ should 
be expanded in a Taylor series around Mq — xi, since 
the diagrams in o-'^"'^ are of 0{p^). If AI^^ is formally 
written as 



'(4) 



(52) 



then an expansion of ct up to 0(ji ) gives 



IV. ANALYTICAL RESULTS FOR D 

Equation H47|l is written in terms of the blocks of , 
but for the explicit calculations we want to have the ex- 
pression in terms of the self-energies Ey directly, since 
this is what we know how to calculate. In the degen- 
erate case we are studying in this paper, there are sev- 
eral equivalent ways to obtain the quantities needed since 
a, p, T, wi and W2 appear in many different elements of 
the full self-energy matrix. We chose to use 

a = Ell — Ei2 , p — 3Ei2 
wi = E44 - E45 , W2 = 3E45 , T = 3Ei4. (48) 

which follows immediately from Eqs. (|3U|I and (|31|1 . Then 
we get 



2?= y(xi-X4 + Eii-4Ei2-E44-2E45 + 6Ei4). (49) 



/(4) 



/(4) 



Mt 



dcT'^^\p',Xr) 



dp^ 



+ 0(/), 



MS 

(53) 



where Mf represents the NLO correction to the lowest or- 
der charged meson mass, which is given by (T*-^''(Mq , Xi)- 
In practice, however, this Taylor expanded term is not 
needed, because cr'''*-' contains no higher powers of p^, 
and thus the derivative term of H53|) is identically zero. 
Collecting the relevant terms up to 0{p^) gives 



Z = 



. + a'^'^ {Mix.) + a'^'^ (A/o^ X.) + {^'^'^ (Mix.)) ' ■ 

(54) 

Next, the form of 2T,ij has to be determined up to 
0{p^). In order to do this, we determine E^ up to 0{p^), 
multiply the expressions for Z and Ey together and col- 
lect the resulting terms up to 0{p^). As for the a terms 
we write the self-energies Ey as a string of terms of pro- 
gressively higher order: 



A. V at NNLO 

For the NNLO calculation of T), we have to determine 
the expression for (|49|l up to 0{p^). This means that all 
quantities appearing in Eq. (|47|l have to be expanded up 
to the appropriate order. 

The first terms are of the form (xi — X-i)^- Here 
Xi — X4 is of 0{p^), so we have to determine Z up 
to 0{p'^). From the definition of Z and a'{M^f^) = 
d(T{p'^,Xi)/dp\2^Ml^, we get 



E,,(m2„x.) = E^f (Mf^xO+E^f (Mf„x.) 



(55) 



For E^^'' the argument must be expanded in a Taylor 
series around Mg = xi; since the diagrams in e|^^ are of 
C'(p^). An expansion of E^^^ up to 0{p^) gives 



^■\mIx.) + 



o(/), 



(56) 



z = 



1 



= l + a'(A/2j+a'(Af2j2 + ... (50) 



The expressions for a'{M'^fJ can be written as a string 
of terms which denote the IPI diagrams of progressively 
higher order. 



-a'^'\M^,,X.) + 



(51) 



where Mf again represents the NLO correction to the 
lowest order charged meson mass, given by ^'^•'(A'/q , Xi)- 
Thus, up to 0{p^), the E^j's are given by 



(Mix. 
''HMlxi. 



aEg)(p2,x. 



dp^ 



0{p'),{57) 
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where the last two terms are of 0{p^) and represent the 
NNLO contributions. For convenience, we will use the 
notation ^'["f for the derivative term, assuming it un- 
derstood that the derivative is as in Eq. H57|) . 

Having this result, we now get the final expression for 
ZEy up to 0{p^) to be 



= sg)(M^,xO + s!f(M^,X.) 



(4), 



(58) 



Putting everything together, the analytical expression 
for V up to 0{p^), or NNLO, is given by 



p(6) 



(Xi - X4) 



-2 E 



+ cr 

(4) 
12 


(4) v'(4) 
^ 11 


+ <7 

(4) , 
12 1 


(4) 
-cr 


+ 

(4)y(4) 
^12 


^11 
+ Sf,' 


(4) 
14 




(4) , 
14 I 


-a 


(4)y(4) 
^14 




(4) 
44 




(4) , 
44 ' 


- cr 


'(4)y(4) 
^44 




(4) 
45 




(4) , 
45 ^ 


-cr 


(4)y(4) 
^45 





(59) 



lowest order meson masses, since these are the funda- 
mental inputs in PQ^PT. The situation is different in 
standard xPT, where the main objective is comparison 
with experiment, in which case the formulas are most 
often rewritten in terms of the physical decay constants 
and masses. We have changed here the general F to Fq 
as is appropriate for the three flavor case. 

The lowest order contribution to V is just the difference 
between the valence sector quark mass and the sea sector 
quark mass, i.e 



V 



(2) 



3^ 



14- 



(61) 

Here one sees explicitly the fact that this is a quantity 
relevant only in the partially quenched theory, since the 
lowest order contribution vanishes when one takes the 
limit X4 ^ Xi of unquenched xPT. We have checked that 
this also holds for the C(p'*) and 0{p^) contributions. 

The combined NLO contribution to V (loops and coun- 
terterms), is 



2?^^^ = ^{-24LlRf^x4-ieLlRf^xi+^8LlR 



d 

liXi 



-48L;«)2 + 32LlRf^xi4 - 10/3A(xi)xi 
+6i(xi4)xi4 - 8/3i(x4)X4 
-l/3S(xi,Xi;0)i?^4Xi}. (62) 



where the arguments of the functions have been sup- 
pressed. 



and is in agreement with Refs. The NNLO contri- 

bution from the 0{p^) counterterms is given by 



B. Expression for T> 

The analytical expression for T> is given to NNLO in 
the form 



I?(4) Vif 



f2 
^0 



(6) 
loop 



p4 



0(/ 



(60) 



where the 0{p*) and 0{p^) contributions have been sep- 
arated. The NNLO contribution I?'^' has been further 
split into the contributions from the chiral loops and from 
the 0{p^) counterterms or LECs. 

It should be noted that we have chosen to give the 
results to the various orders in terms of the three fla- 
vor lowest order decay constant Fq and in terms of the 



(6) _ 



i{ - 96KI,xIRm - 192Kl^xiX4Rfi 
+Kl^[16XiiRf4f - 48x?i?^4] - 96K'2oXiX4Ri4 

-48if2^lX4< - lUK^22XlRf4 

+K^:,[16xi{Rt)' - 48x?i?^4] + 48X2^4Xi«)' 
+48if2^5[x? - X^] + K;,[96xiX4Rii + 144xli??4] 



xl] 



-432if2^X4< + 32/^3^9 [X? 
-192ifIoXi4X4i?M - l92Kl,xi4{Ri 

-288ifl2X4(i?^4)'}- 



(63) 



The NNLO loop contribution is considerably longer. 
Therefore we have made the FORM '2^ output with all 
the analytical expressions available for download from 
the web site 21]. The NNLO loop resuh is 
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Aoop = ^{7ri6LS[49/3xiX4< - 2Xii«M + 3x^4] + 127ri6iiXi< + 7ri6L^[6x?< + 16X4<] 

+7ri6LS[71/6xiX4i?f4 - 8xlRt + V^xIrU + wU21/S2xiX4Rii + 73/6AxIrU + 768L^LgxiX4i?5 
-n52LlLlxlRL + ll?>2LmxMAf ' 768L^Lgxi4X4i?f4 + 576Lf xl^^w - 768igigxiX4iif4 
+384i^i?Xi(^^w)' + LlLl[128xi{Rt^f - SSixMi] + 192Lf x?i?f4 + ^(Xi)^i6[9/4xiX4 + 3/4x?] 
+A(xi)LS[32xiX4 - 80x? - 8x^] + 16A(xi)iIxii?f4 + 40i(xi)Lki^u 
+A(xi)L^ [32x1X4 - 44x? - 8x1] + 160i(xi)ikiX4 + A(xi)LS[-32/3xiX4 + 208/3x?] 
+A(xi)ig[-224xiX4 + 64x?] + 6iA{xiWrXiRf4 - 352/3A(xi)i:gx? + 14/9A(xi)\i 
+7/2^(xi)A(xi4)xi +^(Xi)5(xi,Xi;0) [-10/9x1X4 + 10/3X?] - 2^(xi)5(xi4,Xi4;0)xf4 
+2/9A(xi)C'(xi,Xi,Xi;0)x?^u + ^(xi4)7ri6[-3/4xii?f4 - 3X4] + A(xi4)i:o[144xiXi4 - 24xi4X4] 
+i(Xi4)i^3[144xiXi4 - 60X14X4] - 288l(xi4)i4Xi4X4 - 144A(xi4)i5XiXi4 + 288A(xi4)i^6Xi4X4 



^288^(Xi4)i;Xi4i?u + A(xi4)i8[192xiXi4 + 86x14X4] - 27/2A(xi4)'xi4 + 8A(xi4)A(x4)X4 
-2A(xi4)S(xi, Xi; 0)xiXi4 - 64A(x4)iSxiX4 + 128 A{x4)L{x4Rti + 32A(x4)ik4i?u - 64A(x4)ikiX4 
+i(X4)iI[- 128x1X4 + 256x1] + 256/3A(x4)iki4X4 + i(x4)i^[128xiX4 - 256x1] + 256i(x4)i^X4i?w 
-512/3^(X4)ik4 + 4/9^(x4)'x4 + 8/9^(x4)5(xi, Xi; 0)xiX4 + 8/9^(x4)S(X4, X4; 0)xl 
+A(xi, e)^i6 [7/36x1X4 + ll/18x? + 1/4x4] + ^(Xi4, e)7ri6[-19/2xiX4 - 5/2x? + 9/2x1] 

+58/9A(x4, e)^i6X4 + 5(xi, xi; o)LS[8xiX4< " ^^xIrU + 8S(xi, Xi; OWixlR^ 

+20B(xi,Xi;0)£k?< + S(xi,Xi;0)LS[8xiX4i?M - 16x?<] +^(Xi,Xi;0)i^[- 24x1x1 + 104x?X4] 

+B(Xi, Xi; 0)L^[16/3xiX4 - 80/3x?X4 + 48x?] + S(xi, Xi; 0)Lg[32xi«)' - 160x?X4] 

+48S(xi,Xi;0)i?Xi(i«?4)' + B(xi,Xi;0)Lg[16/3xixl + 64/3x?X4 - 80x?] 

+B(xi, Xi; 0)2[-l/18xiX4< + 23/18x?i??4] + 2/9B(xi, Xi; 0)C(xi, Xi, Xi; 0)x?(i?f4)' 

-144B(xi4, X14; 0)L^x?4X4 - 48B(xi4, Xw; 0)L^x?4 + 288B(xi4, Xi4; 0)i6xLx4 + 96B(xi4, X14; 0)iSx?4 

+64B(x4, X4; 0)L^X4 + 64/3B(x4, X4; 0)Llxl " 128B(x4, X4; 0)Lgx4 - 128/3B(x4, X4; 0)L^X4 

+B(X1, Xi; 0, e)^i6[-l/4xiX4i?f4 " 2/3x?<] + 16C(xi, Xi, Xi; 0)Lk?X4ii?4 

+16/3C(xi , XI , XI ; 0W5xlRf4 - 32C'(xi , XI , XI ; 0W6xlx4Rf4 - 32/3C'(xi , XI , XI ; 0)LlxlRf4 

+H^{1, XI, XI, XI, XI ) [-4/27x1X4 + 44/27x?] + 5/9ff^'(l, Xi, Xi, Xi, xOx?^?^ 

+fl'^(l, XI, Xi4, Xi4, Xi)[-10x? + l/4X4i?u] + H^' (1, XI, Xi4, Xi4, Xi)[l/4xiX4i2u - X^M 

-2H^{1, xi4, Xi4, X4, Xi)X4 + 2-ff^' (1, Xi4, Xi4, X4, Xi)XiX4^f4 - 40/27i?^(l, X4, X4, X4, Xi)xl 

+5/3F^(2, XI, XI, XI, XOX^^-M + 4/9i/^'(2, Xi, Xi, Xi, Xi)x?(^«f4)' 

+H^{2, xi, Xi4, Xi4, Xi) [1/4X1X4 i?f4 - ^Ma] - 5/4if^'(2, xi. Xi4. Xi4, Xi)x?(i?u)' 

+4/97?^(5, XI, Xi, Xi, Xi)x?«)' + 2/9if^'(5, Xi, Xi, Xi, Xi)x?«)' 

+2/27iJ^(8, XI, XI, XI, Xi)x?(^w)' + H[{1, xi, Xi4, Xi4, Xi)[2xiX4 + 18x?] 

+16iJf (1, X4, Xi4, Xi4, Xi)XiX4 + 2H[ (2, xi, Xi4, Xi4, Xi)Xi(-Ri4)^ + 2iJf (2, Xi, Xi4, Xm, Xi)Xi(^i4)^ 
+F2''i(l,Xi,Xi4,Xi4,Xi)[- 3/2x1X4 - 6xl]+3/4Hi^ (1 , Xl , Xl4, Xl4, Xl)Xl-Rl4 

+27/2Ffi(l, xi4, Xl, Xi4, Xi)x? + -fffi(l, X4, Xi4, Xi4, Xi)[-12xiX4 + 6x?] + 6Ff/(l, X4, Xi4, Xi4, Xi)x?-Ru 
+Ffi(2,xi,Xi4,Xi4,Xi)[3/2xiX4Jif4 - 3/4x^4] - 3/4Ff;(2, xi, Xi4, Xi4, Xi)x?(^«f4)'}- (64) 



V. NUMERICAL RESULTS AND DISCUSSION 

A. Checks on the Calculation 

As for the previous articles on NNLO calculations in 
PQxPT, the calculations are rather lengthy, and require 



careful checking procedures to avoid the possibility for 
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mistakes showing up in the end result. The complete cal- 
culations have therefore been performed in two complete 
independent versions, which have been checked against 
each other for agreement. The numerical results have 
been treated likewise, and the programs were written 
in different programming languages (Fortran and C++). 
Furthermore, the divergence structure agrees with gen- 
eral calculation of J^] and the cancellation of nonlocal 
divergences happens as required. The end results are ex- 
pected to vanish for X4 ~^ Xii since the quantity V is 
a quantity relevant for the PQ theories only. This has 
been shown both for the O(p^) and 0{p^) analytical for- 
mulae and checked numerically as well. The latter can 
also be seen from the fact that the contour plots below 
do not show a peaked behavior near the diagonal where 
X4 = Xi- 

We have checked explicitly that the self-energy when 
calculated for the different channels has the structure of 
Eqs. l|5n)l and as well as the relation of a to the 
charged case of Ref . . 




0.05 0.1 0.15 0.2 0.25 0.3 
Xi [GeV^] 



B. The Numerics 

All plots in this section are given in terms of the rela- 
tive shifts to the lowest order contribution 'D^^\ That is, 
we plot 



where we keep the contributions in Eq. H60I) up to and 
including 0{p'^) when we refer to the NLO case and up 
to and including 0{p^) when we refer to the NNLO case. 

In the long run, the input parameters should be de- 
termined by fits of the PQxPT formulas to lattice QCD 
data. These fits are not available at present, although 
suitable simulation results should become available in 
the near future. So we present here plots with the same 
choices of input parameters as in JT^ JJ^, . This is 
mainly from the NNLO fit to experimental data, referred 
to as "Fit 10" , which is presented in Ref. . That fit has 
Fq = 87.7 MeV and a renormalization scale of /i = 770 
MeV. The NNLO LECs Kl and the NLO LECs L^, 
and Lg were not determined in that fit and they have 
been set to zero for simplicity. It should be noted that 
Lq cannot be determined from experimental data, since 
it is a distinguishable quantity only in the PQ theory. 
Some recent results on L4 and Lg have been obtained in 
Ref. "251, but they have nevertheless been set to zero in 
the plots in this paper, since the present numerics are 
mainly intended for illustrative purposes. 

For the degenerate mass-case we are working with in 
this article, we have only two input parameters for the 
quark masses, namely xi ^'^d Xi- This allows for ex- 
haustive contour plots of the whole parameter region of 
interest. A contour plot of the NLO relative correction 
is given in Fig. |21 and in Fig. |31 the same parameter re- 
gion is plotted for the sum of the NLO and the NNLO 



FIG. 2: The NLO relative shift of the lowest order contribu- 
tion to the double pole coefficient AT), as a function of 
the valence and sea-quark masses xi a-nd X4- The difference 
between two successive contour lines in the plots is 0.10 and 
the values chosen for the LECs correspond to "Fit 10", as 
discussed in the text. 



shifts. Let us remind the reader here that the quantities 
Xi = 2Bomi correspond to lowest order meson masses 
squared. For the NNLO case, diagrams appear which 
can have on-shell intermediate states for xi > 9X4- For 
that part of the plot we have simply kept the real part 
only of the relevant sunset integral. The error due to this 
should be a small effect, similar in size to the width com- 
pared to the mass of the eta for realistic quark masses. 

The plots show a reasonable convergence behavior for 
when both XiiX4 ^ 0- Furthermore, as mentioned 
above, since we have plotted the relative correction with 
respect to the lowest order contribution I?'^' = Xi ^ X4i 
the plots would have a peaked limiting behavior near the 
diagonal line xi = X4 if NLO or NNLO corrections 
were nonzero for xi — Xi- That this is not the case is 
clear from the plots, and in agreement with the expected 
result. Furthermore, the general shape of the contour 
over the plotted parameter space is rather different for 
the NLO and the NNLO results, in particular, the ©(p*^) 
contributions are rather sizable. This is an expected re- 
sult as well, since previous calculations of the masses and 
decay constants to NNLO show exactly this behavior. 
In Ref. it was shown that the type of behavior is 
strongly dependent on the values of the largely unknown 
NNLO LECs and an example of changes in some of the 
less well determined LECs was given that improves the 
convergence behavior dramatically. 

However, since a major point of this article is the de- 
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AD, Fit 10, 




Xi [GeV^] 



FIG. 3: The combined NLO and NNLO relative shifts of the 
lowest order contribution to the double pole coefficient D, 
AD, as a function of the valence and sea-quark masses xi and 
X4. The difference between two successive contour lines in the 
plots is 0.10 and the values chosen for the LECs correspond 
to "Fit 10", as discussed in the text. 



0.2 




-0.7 I ' ' ' ' ' 1 

0.05 0.1 0.15 0.2 0.25 0.3 

Xi [Gev2] 

FIG. 4: The NLO relative shift of the lowest order contribu- 
tion to the double pole coefficient V, AC, plotted for 6 — 60° 
in the xi - Xi plane, for 10^ = {0, -1, -2, -3, -4, -5, -6} 
(top to bottom curve). The values chosen for the other LECs 
correspond to "Fit 10" as discussed in the text. 



pendence on the LEG for D, we have, in addition 
to the contour plots, chosen a parameter hne in the two- 
dimensional mass-space, where results for different values 
of the LECs have been plotted. This parameter line is 
defined by 

X4 = tan60°xi, (66) 

This choice is motivated by the fact that the sea quark 
mass is usually taken to be larger than the valence quark 
mass in lattice simulations. 

In Fig. ^ the NLO relative correction has been plotted 
along this parameter line, as a function of xi- The plot 
is for several different values of the LEG Ly, ranging be- 
tween zero and —0.0006. All other parameters are set to 
the values of "Fit 10". Convergence as xi 0, for all 
values of Uj, is evident, in agreement with the contour 
plot. The dependence on is exactly linear as expected 
from 9] and Eq. (|^ . 

In Fig. 13 the NNLO relative correction is plotted us- 
ing the same values for the LEGs, namely "fit 10", ex- 
cept for Uj. As a comparison. Fig. shows the NNLO 
relative correction with all LECs, except Ly, set to zero. 
As anticipated from the contour plot, the contribution is 
sizable in both cases, but with a somewhat lower overall 
amplitude when the LEGs are set to zero. Also for the 
NNLO results, we have an exact linear dependence on 
as can be seen from Eq. H64|) . As already discussed, the 
large amplitude of the NNLO corrections is not expected 



1.4 



Q 




-0.2 I ' ' ' ' ' 1 

0.05 0.1 0.15 0.2 0.25 0.3 

Xi [GeV^] 

FIG. 5: The NNLO relative shift to the lowest order contribu- 
tion to the double pole coefficient D, AD, plotted for 9 — 60° 
in the xi - X4 plane, for 10* = {0, -1, -2, -3, -4, -5, -6} 
(top to bottom curve). The values chosen for the other LECs 
correspond to "Fit 10". 
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-0.2 I ' ' ' ' ' 1 
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FIG. 6: The NNLO relative shift to the lowest order contribu- 
tion to the double pole coefficient AD, plotted for 9 = 60° 
in the xi - Xi plane, for 10* L? = {0, -1, -2, -3, -4, -5, -6} 
(top to bottom curve). The values chosen for the other LECs 
are all set to zero. 



to be a concern for the convergence of the theory, since 
many of the LECs are as of yet poorly determined. As 
an example to show this as well as the subtraction scale 
dependence, we have shown in Fig.[7|the NNLO relative 
correction as a function of iy and all other LECs set to 
zero at a subtraction scale /i = 0.6 GeV. 



VI. FITTING STRATEGIES 
A. LECs relevant for the eta mass 

The eta mass to NNLO in xPT was calculated in [T9l | 
and to NLO in [sj. We present here the corrections in 
the form with the dependence written in terms of lowest 
order masses and decay constants as calculated in [T^ 
and obtainable from The NLO part depending on 
the p-* LECs is 

m2(4) ^ TOm,(-40/3X^-32/9L^-|-80/3i:^-64/3i;) 
+m^{-lQ/ZLl - 8/9LI + 32/3Lg + 32/3i; 
+IG/3LI) + ml{~m/3Ll - 32/9i^ + 32/3Lg 
+32/3^ + 32/3LI) . (67) 

The occur in the three combinations 2Lg— L4, 2_Lg— 
and LI + 6L^. 
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FIG. 7: The NNLO relative shift to the lowest order contribu- 
tion to the double pole coefficient V, AC, plotted for 6 — 60° 
in the xi - Xi plane, for 10* = {0, -1, -2, -3, -4, -5, -6} 
(top to bottom curve). The values chosen for the other LECs 
are all set to zero. The subtraction scale is /i = 0.6 GeV. 



The NNLO part depending on the NNLO LECs is 



to2(6) _ TOm2(-128/9C[2 - I28/3Q3 - 32/9C[4 

-64/3C[5 - 16/3C[g - 32/9C[7 + 32/3C{^ 

+48^2^0 + 144^2^ + 128/3^3^2 " 64/3CJ3) 
+rh'^ms{-U/9Cl^ - 32C[3 - 32/9C[4 
-16C[5 - 64/3Cr6 - 32/9Cr7 + 32Cl„ + 192C2''i 
+32/3CJ2 - 64/3CJ3) + m3(-32/27q2 
-64/9C[3 - 16/9C[4 - 32/9C[5 - 32/3C[6 
-16/9C[7 - 32/9C[8 + 16C[9 + 32C^o + 64C^i 
+32/3CJi + 64/3C3^2 + 64/3CJ3) 
-l-m3(-256/27C[2 - 128/9C[3 - 64/9C[4 
-64/9C[5 - 32/3C[g - 64/9C[7 - 64/9C[8 
-H32C[g -f 32C2^o -I- 32C^i + 64/3CJi -f 64/3C3^2 
+64/3CJ3) . (68) 



The CI contributions there can be rewritten in terms 
of the three flavor PQxPT NNLO LECs Kl using the 
Cayley-Hamilton relations of [ItII . The arguments for this 
were presented in detail in |l3l | . The resulting formula in 
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terms of tlie combinations of Kl defined below is 

m2(6) ^ rhml{~i2/9X2 - 128/3X3 - 48^4 + 128/3^5 
-64/9^6) + m2m,(-32/9X2 - 32^3 - 64^4 
+128/3X5 - 64/9X6 + 64/9X7) + m3(-32/3Xi 
+80/9X2 - 64/9X3 - 64/3X4 + 32/3X5 
+64/9X6 - 128/27X7) + m3(-64/3Xi 
+128/9X2 - 128/9X3 - 32/3X4 + 64/9X6 
-64/27X7) . (69) 



B. Obtaining LECs from V: fitting strategies 

In Ref. \\^ we discussed at some length, possible fit- 
ting strategies for determining the various LECs at Olp^) 
and 0{p^). In the analytical expressions for V some ad- 
ditional LECs are present, which motivates further elab- 
oration on these fitting strategies. At 0{p'^) the only new 
LEC is Lj and at 0{p^) we now have i<r24j ^41 ^-'^d 
as well. It turns out that together with the expressions 
for T) it is possible to determine all the combinations of 
LECs which appear in the NLO+NNLO expressions for 
the masses including the eta mass. We will now give a 
short outline of how one might proceed but we will con- 
centrate on the masses only. 

From the charged masses we have the dependence 
J(f ' ~ 3x m - Ll) + X13 (2L^ - LI) . (70) 

In these expressions, x denotes the average over the sea- 
sector quarks, and the superscripts 23 indicate having 
2 nondegenerate valence quark masses and 3 nondegen- 
erate sea quark masses. We refer to the discussion in 
Ref. ^3 for details on which mass degeneracies which 
are needed for the determination of the LECs, and just 
refer to the needed expressions here. 

The NLO result for T) is proportional to 

~ 3x4(2L^ - LI) + 2x14(22.;; - Ll) - R^Ll + 6L;) 

(71) 

from which an overall factor of Rf^ has been removed. 
The resemblance to the mass expressions is striking, even 
though the coefficients are not quite the same. However, 
it is clear that all combinations needed to determine the 
eta mass at NLO can now be obtained. Since we can 
obtain L5 from the decay constant we can also get at L7 
separately. This was already shown to NLO in jg, ■ 

Let us now do the same thing for the NNLO LECs. The 
situation is somewhat more involved but not different in 
principle. The full discussion can be found in [T^. we 
only present a part of it here. The mass at NNLO for 



the (2+2) case there depended on the NNLO LECs as 

= XiX3X4(-32X3) + XiX3X6(-16X3) 

+XiX3(-12^i + 8X2) + xiX4X6(-32X4 + 32X5) 
+XiX4(-32X4 + I6X5) + xiX6(-8^4) 
+X?X3(-12Xi + 8X2) + x?X4(-16X3) 
+X?X6(-8X3)+x?(-4Xi) 

+X3X4X632(-X4 + X5) + X3X4(-32X4 + I6X5) 
+X3xi(-8X4) + X3X4(-16X3) + X3X6(-8X3) 
+X3(-4^i)- (72) 

The combinations Xi are defined as 



Xi 


— -"^39 






X2 


= Kl, 


K23 ^ 31^25 




^3 


= Kl,+ 






X4 


= ^22 + 


K21 — — 32^27 




^5 


= ^21 - 


-"■26 




^6 


= ^23 + 


2KIj - i/2K^2i + 3/2iir25 


+ 3Kl^ 


X7 


= K'u- 


3/2i^2-4 - iKl2, - 3Kl^ - 


9KI2. (73) 



From Eq. (|72() it is clear that the combinations Xi to X5 
can be determined from the charged masses at NNLO 
by varying the input quark masses in the lattice QCD 
simulations. The above is a rephrasing of the discussion 
of fisl. 

Let us now rewrite the NNLO LECs contribution to V 
in the same fashion. Eq. 163|) can be rewritten as 

T^if ~ XiX4(-2Xi+X2 -12X3 + 2X7) +x?(-2Xi 
-2X6) + xI(-2Xi + 2X2 - 9X4 + 6X5 + 2X6 
-2X7) , (74) 

It is clear that from V wc can thus determine the remain- 
ing LECs needed to get at the eta mass at NNLO. We 
can easily get at the missing combinations Xe and X7. 

VII. CONCLUSIONS 

In this paper, we have calculated the quantity V to 
NNLO. This quantity was first proposed and calculated 
to NLO by Sharpe and Shoresh I3 , as a way to obtain 
a better fit of the LEC L7, since it appears already at 
NLO for this quantity. This then also allowed to get the 
eta mass at NLO in xPT. As discussed in 's'l, I? is a 
well defined Lattice quantity in terms of the correlation 
function Ro{t). 

We have presented the full analytical results for V to 
NNLO, together with a numerical analysis indicating the 
typical size of the corrections for various choices of the 
LECs. In particular, we demonstrate the sensitive de- 
pendence on the LEC L7 which these expressions ex- 
hibit, both at NLO and at NNLO. It should also be noted 
that, in addition to showing dependence on the LEC Lj 
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already at NLO, the NNLO analytical expressions pre- 
sented here also depend on the LECs K24, and K\2i 
none of which were present in the NNLO expressions for 
the charged masses or decay constants. 

We have extended the discussion about possible fitting 
strategies in Ref.0| to include the new LECs which ap- 
pear in the NNLO results for T). The conclusion is that 
the NLO and NNLO expressions for the charged masses 
and decay constants, together with the expressions for 
T) provided in this article, in fact give enough informa- 
tion to determine all of the combinations of LECs which 
appear in these calculations, including those for the eta 
mass at NNLO. 

In addition we have also shown how to resum explicitly 
the full propagator from the lowest order propagator and 
self-energies also in the case for the neutral propagator 
in PQxPT. 
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APPENDIX A: SYMMETRY CONSTRAINTS ON 
THE PROPAGATOR G^jki 

In Ref. 1^, App. E, a number of constraints on the 
form of the meson propagators in PQQCD were de- 
rived. The Lagrangian of QCD played no role, since the 
only properties of the quark bilinears which were used 
were the transformation properties under vector symme- 
tries, whose properties are shared by the meson fields of 
PQyPT. It was then argued that the derivation presented 
in Q could easily be extended to the case with more than 
one flavor of valence quarks and nondegenerate masses. 
We will now show this generalization to the general case 
with rinvai valence and ghost quarks and risea sea quarks. 
The masses are fully general except that the a-th ghost 
quark has the same mass as the a-th valence quark mass 
for all a — 1 , . . . , Jinvai • Three types of symmetries are 
needed, all fairly trivial extensions of the ones in 

We suppress the time ordering and often the space- 
time dependence in the remainder of this appendix to 
make the symmetry structures more visible. 



1. Phase rotations and G"'^ = G" 
We write the quark flavors as a vector 

q = (gi, . . . , g2n,„,,+«„,ea)'^> (Al) 



where indices 1 to rinvai denote the valence quarks. In- 
dices rinvai -I- 1 to rinvai + "-nsca dcnotc the sca quarks 
and indices Unvai + "•nsca + 1 to 2nnvai + "-nsca denote the 
bosonic "ghost" quarks. 

The direct generalization of the subgroup V of the vec- 
tor transformations in Ref. j^, is 

V = diag (exp 6*1 , . . . , exp 6'2n„™i+n„3ca ) ■ (A2) 

Under a phase rotation of only the flavor to, Gijki trans- 
forms as 

Gijki Gijki exp {6„i{Sim — Sj„i + 5km ~ Sim)) ■ (A3) 

The invariance of Gijki implies that the indices must be 
paired up. Thus, non-vanishing elements of G have the 
form Giijj or G^jji . 

A second relation which can be simply derived is that 
G" is a symmetric matrix. Consider 

($y(x)$fc,(o)) = {^.,{~x)^kim 

= ($y(0)$fci(x)) 

= (-)(<&,, (^)<i>.,(0)), (A4) 

which make use of rotation invariance and translation 
invariance. The minus sign in the last step is only needed 
when both $jj and are fermionic fields. This equation 
implies that Gijki = {—)Gkiij- In particular, this implies 
that G" defined in \i'A\ is a symmetric matrix. 
At this level G" has the structure 

(r -f s F \ 
t u w \ . (A5) 
V vo^ r ~ s I 

r, s,u are symmetric matrices. The matrices are ordered 
in the valence, sea and bosonic sectors. 

2. SI7(1|1) transformations 

The third type of symmetries is SU transforma- 
tions involving the TO-th quark from the valence sector qy 
and the m-th quark from the bosonic sector qi,. There are 
rival of such SU {1\1) symmetries. Below we will use the 
label niv and tob to denote these quarks respectively. 

This type of transformations is of the form 

U^l^l J) eSU(l|l)„. (A6) 

where a, d are commuting numbers, b, c are anticommut- 
ing numbers and U satisfies UW = I2x2- Here we think 
of it as acting on the vector 

Q = {qmv , IniB y ■ (A7) 

The index to in (jA6|) indicates which S't/(1|1) we are 
referring to. 
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This type of symmetry transformations can, together 
with the previous symmetries, be used to set further con- 
straints on the block structure of Gtjki- In order to see 
this, we first form the following classes of 2-indexed ob- 
jects out of elements of G: 



^ik 



E 



(A8) 



where ns indicates the n-th index in the sea sector, 
and the matrix indices i, k take both take values from 
{my, ms} so as to form a supermatrix (or graded) struc- 
ture. The symbol ej gives the correct sign of the graded 
symmetry, i.e. ej = 1 for j = l,...,nvai + ^sca and 



Ej = —1 for the remaining ones. The objects in( IA8|) 
are the direct generalizations of the objects in 9] to the 
general flavor case. 

In addition to these, we will also need the two classes 
of objects defined by 



n° — c 

6 
ik 



/o6 _ (-1 
^ik ~ '^nsni 



ik{x) , 

ik{x) ■ 



(A9) 



Here, the value chosen for ny must be different from the 
value chosen for my and ub different from m^. The 
subscripts V and B indicate that ny is a valence index 
and ub & ghost index. All the objects defined above 
transform imder SU{l\l)m as 



o: 



ik 



/ , ^lo'^op^pk 



(AlO) 



o.p—my ,m£ 



Since SU{l\l)ra is a symmetry, the expressions are in- 
variant under these transformations. Writing out these 
transformations explicitly one arrives, after some calcu- 
lations, at a number of constraint equations, which must 
be satisfied, for the symmetry to hold. Solving these 
equations, the constraints following from (jA8|l are 



G, 

a 

G. 



mv Tny rny my 



mvmYmBmE 



niB'iT^vfnY'n^E 



'isnsmvmv 





Gfn BrriBiTtvmv 

-G. 
G 



'niYfnBfnBfnY 



^mvrnvnsns 



— Gnsn 



snsmBiriB 



G 



mvnsnsmv 



GfrLBf^Bf^s'TLs — ^mn 
GrLsmvmvns — GfyiBn-snsiriB 



— —Gns 



srriBfnBns 7 



(All) 



where as defined above, my is the m-th valence quark 
index m^ is the m-th ghost quark index and ns is the 
n-th sea quark index. What we have shown so far is that 
the matrices t and w in (|A5|I are identical and that the 
diagonal parts of r and v are the same. 



The previous was a simple extension of what was done 
in ,9] . The really new part is that for the general flavor 
case, we also have a distinction between diagonal ele- 
ments and non-diagonal elements in the various matrices 
of G". The invariant objects in (|A9|) provides informa- 
tion about the matrix structure of these. Writing out 
the transformation (|A10|I for these objects yields the ad- 
ditional constraint equations 



Gnvn 



G 



vnvirLvmv 



GfLYTL 



YnvmBrriB 



riBnBiTiYfnv 



G 



(A12) 



for ny 7^ my, ns ^ tub- These equations imply that the 
off-diagonal part of s vanishes and that the off-diagonal 
parts of V and r are the same. 

To summarize, G" has the structure 



G" 




(A13) 



r,s,u are symmetric matrices and s is diagonal. 

Here we can see that since t and r connect different 
sectors they are disconnected if we think in terms of ex- 
ternal quark lines. That s is from connected lines and 
related to charged propagators is a little more involved 
to show. First, we use the fourth line in Eq. (|Alip where 
we have an obviously charged object in the propagator. 
We now set a second valence/ghost quark to the same 
mass, i.e. rrimv — frinv ■ We can then use the permuta- 
tion symmetry to turn the my indices into ny . The final 
step is to use the SU{l\l)m symmetry to rotate the m^ 
index back to my. 

^nun — GmBnivmv^B GmBnvnv'tnB Gm^n^nv^^^v? 

(A14) 

which is the desired relation between s and the charged 
meson propagator. 



3. Equal mass limits 

When masses are equal, we have more symmetries. 

For all valence quark masses equal we can interchange 
two indices my and ny everywhere and G" should re- 
main invariant. This implies that we can interchange 
columns in t and v of Eq. ljA5p . It also implies that r + s 
must be invariant under the simultaneous change of a 
pair of columns and rows. 

For all bosonic quark masses equal we can interchange 
two indices tub and ub everywhere and G" should re- 
main invariant. This implies that we can interchange 
rows in v and w'^ of Eq. ljA5(l . It also implies that r — s 
must be invariant under the simultaneous change of a 
pair of columns and rows. 

For all sea quark masses equal we can interchange two 
indices mg and ns- This means that t and w are invariant 
under the interchange of rows and that u is invariant 
under the simultaneous interchange of rows and columns. 
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Putting the above together with t = w and r = v 
derived earher shows that r and t are matrices with all 
elements equal, s was diagonal so all its diagonal ele- 
ments must be equal. For u it implies that all diagonal 
elements are equal and that all off-diagonal elements are 
equal. The matrix then has precisely the structure of 
Eq. 



APPENDIX B: THE RESUMMATION IN THE 
GENERAL CASE 

In this appendix we show how to do the resummation 
in the general case for an arbitrary number of valence 
quarks rival and sea quarks ngca- 

The general lowest order neutral propagator matrix 
(|16|l can be written in the form 



a 
G° = i\ f3 \-iS 
-a 




(Bl) 



The notation here is similar but not the same as in the 
main text. The matrices a and f3 are diagonal matrices 
with propagators. 

a = diag(ai,...,a„^^j) , 
/? = diag (/?!,..., /?„^^ J , 
5^ = (ai, ■ • ■ ,anvai) I 

= (/?!,..., , 

"i = 1/ {P^ - Xi) , 
5 = 



1 



(B2) 



The the various a symbols are the valence and ghost 
sectors and [3 the sea quark sector. In the equal mass 
limit Eq. (jBTJ reduces to Eq. ((T7jl . 

The self-energy has the structure derived in App. ^ 
where the submatrices now have a more general structure 
than in the main text. 



R + S -R 
E = I T W -T 

-R -T^ R-S 



(B3) 



i? is a symmetric rival by rival matrix. is a symmetric 
JT-soa by risca matrix. S' is a diagonal rival by rival matrix. T 
IS a rigQa by rival matrix. R and T have only disconnected, 
S only connected parts when thinking in terms of the 
external quark lines. W has contributions from both. 
The full propagator in the neutral sector is the sum 



G= J2 (G°H)S)"G°. 



(B4) 



n— 0,00 



Not all parts in S will appear in the final resummation to 
all powers. We thus need to extract the parts that appear 



to all orders similar to what was done in the main text 
by splitting into I and K. We do it here by splitting 
(G*'(— i)5]) into four parts with different properties. 



G"(-i)S ^A + B + C + D, 



(B5) 



with 



aS 
A = \ (3W 
aS 

aR -aR 
B = \ 
aR ~aR 

aT^ 
C ^ \ (3T -(3T 
aT^ 

-St 

D = I -Si3 \{a'S + 0''T 
-Sa 



(B6) 



We define four general types of structures 



II 



III 



IV 




U V — U 



(B7) 



The matrices U, V, Y, Z and column vectors u, v, x, y 
stand here for a generic form. These various types satisfy 
matrix multiplication properties. 

• A product of a type I matrix with another one of 
type I,II,III,IV is of type I,II,III,IV respectively. 
This is true for multiplication by type I from front 
or back. 

• A product of a type II matrix with another one of 
type II, III, IV is zero. This is true for multiplication 
by type II from front or back. 

• A product of two type III matrices is of type II. 

• A product of type IV with type I, III, IV is of type 
IV. 

The first property allows to treat the A part similarly 
to the I in the main text. We resum it separately and 
define 



A: 



(B8) 



n— O.oo 
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Due to the diagonal structure of a and /3 this resumma- 
tion can be done with 



A 












• 




• 














1-aS 









1 



(B9) 



l-aS 



The second property means that B can only appear 
once, and together with the third means that C multi- 
plied by itself or products of A appears at most twice. 
Putting then all strings of A only together in A the sum 
can be rewritten as 

J2 (-«G"S)" = X! ((^ + ^CA)D)" {A + Ac A) 

n— 0,00 n — l,oo 

+A + Aba + Ac A + Ac Ac A . 

(BIO) 

The remaining summation in IjBlOjl can be done by notic- 
ing that 

{{A + ACA)Df = {{A + ACA)D) (-J/3^M^/3/3) . 

(Bll) 

The quantity in the last bracket is a scalar and can be 
taken out of the matrix products in the partial sum. The 
sum then can be done by a geometric series in only the 
quantity in the last set of brackets. 

The remaining part of the calculation is to perform all 
of the matrix multiplications. Some properties that are 
used heavily to bring the result into its final form are 



1 



1-PW 



P = P' 



1 



1 - wp 



1 



1- f3W 



1 



(B12) 



There are similar properties for the a and S combina- 
tions. 

The result can be written in a form reminiscent of the 
lowest order (|Blll . 



G = ^ I t w t \ -iS \ b \ (5.^ b'^ d.^) , 

(B13) 

where the matrices r, 3,1,11! have the same structure as 
the matrices R,S,T,U in (|B3|I . These are given in terms 
of the self-energy and lowest order propagator parts as 


















{a- 


r — s 1 









s 
w 
r 
t 
a 
b 

S 



aa , 
P/3, 

aaRaa + aal'^ (ifiTaa , 
aa + aT'^PP , 
6 

1 + sp^^wpp ■ 



(B14) 



Eqs. ljB13|l and (|B14|I are the main results of this ap- 
pendix. Note that all lowest order poles have been shifted 
to the resummed poles everywhere where they appear. 
We have also checked that the results of this appendix 
agree with those in the main text by taking the various 
equal mass limits. 
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